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Abstract 

We study the infra-red limit of non-abelian Chern-Simons gauge theory perturbed 
by a non-topological, albeit gauge invariant, mass term. It is shown that, in this limit, 
we may construct an infinite class of integrable quantum mechanical models which, 
for the case of SU(2) group, are labelled by the angular momentum eigenvalue. The 
first non-trivial example in this class is obtained for the triplet representation and it 
physically describes the gauge invariant coupling of a non-abelian Chern-Simons particle 
with a particle moving on S 3 - the SU(2) group manifold. In addition to this, the model 
has a fascinating resemblance to the Landau problem and may be regarded as a non- 
abelian and a non-linear generalisation of the same defined on the three-sphere with the 
uniform magnetic field replaced by an angular momentum field. We explicitly solve for 
some eigenstates of this model in a closed form in terms of some generalised orthogonal 
polynomials. In the process, we unravel some startling connections with Anderson's 
chain models which are important in the study of disordered systems in condensed 
matter physics. We also sketch a method which allows us, in principle, to find the 
energy eigenvalues corresponding to the above eigenstates of the theory if the Lyapunov 
exponents of the transfer matrix of the infinite chain model involved are known. 
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Gauge theories with Chern-Simons terms [1] have spanned a wide range of interests 
in the past decade or so. The impressive array of topics making use of ideas related to 
Chern-Simons gauge theories extends, on the one hand, from the purely mathematical - 
pertaining to issues of the topology of three manifolds [2] - to the phenomenological, on 
the other, as in models of quantum Hall effect [3] which are testable in the laboratory. 
From a more formal point of view, Chern-Simons gauge theories have shed light on 
aspects of anyonic spin and statistics, and conformally invariant quantum field theories 
in two space dimensions [4]. Driven by a desire to understand some aspects of such 
theories in the simpler setting of 0+1 dimensions, models of abelian Chern-Simons 
quantum mechanics have also been constructed, canonically quantised, and solved [5, 
6] . These models are described by Lagrangians which are quantum mechanical analogues 
of the so-called self-dual models in 2+1 dimensions with and without a Maxwell term 
[7]. The equation of motion of one of these models is given by the famous Lorentz 
equation for a charged particle moving in an external electromagnetic field. As is well 
known, the spectrum of this model is described by Landau levels [8] . The other model, 
obtained by tuning a certain dimensionful coupling to zero, corresponds to projecting 
on to the lowest Landau level. As such, both the models are of immediate relevance in 
quantum Hall effect. 

The conventional way of thinking of the Chern-Simons term is to regard it as a gauge 
invariant mass term in the lagrangian density for a gauge theory in 2+1 dimensions 
with the kinetic piece being given, as usual, by the Maxwell or the Yang-Mills term [1] . 
Contrary to this, one may also construct theories in which the Chern-Simons term plays 
the role of the kinetic term while a Proca term takes the place of a mass term [7]. This 
is the abelian self-dual model mentioned above. As described in [6], the infra-red limit 
of such a theory yields one of the Chern-Simons quantum mechanics models discussed 
in [5]. In the present paper we consider a non-abelian extension of this model and show 

2 



that we may construct an infinite class of quantum integrable models which, for the 
case of SU(2) group, are labelled by the eigenvalue of the angular momentum operator. 
The first non-trivial model in this class will be discussed in some detail and explicit 
expressions will be worked out for some of its eigenstates. We will also sketch a method 
which allows us to calculate, in principle, the corresponding energy eigenvalues. 

To begin with, let us recall some basic notions and establish notations. The non- 
abelian Chern-Simons gauge theory is defined by the lagrangian density 

Ccs = ^ X Tr(A^d v A x + i\a^A v A x ) (1) 

where A^ = A^T a , T a being traceless, hermitian matrices in the fundamental represen- 
tation of the gauge group SU(N). The T a satisfy the following algebra: 

[T a ,T & ]_ = if c ab T c , [T a ,T b ]+ = d c ab T c + ±5 ab , Tr(T a T b ) = 6 -f 

The metric h^ v = diag(— 1, 1, 1), the completely antisymmetric Levi-Civita tensor is 
such that e°* J = e* J = — with e 012 = 1 and we sum over repeated indices without 
comment. ^From simple dimensional considerations it follows that A M has dimensions of 
mass and, consequently, the coupling k is dimensionless. Under a gauge transformation, 

A^x) - A%(x) = U- 1 (x)A^x)U(x) - iU~ 1 (x)d l ,U(x) (2) 

where U(x) G SU(N), the action transforms as follows: 

S[A] -> S[A U ] = S[A] + 2>Kkn (3a) 

where 

n = ^2 / d3x e' M/X Tr(U- 1 d fl UU- 1 d„UU- 1 dxU) (36) 
is the winding number of the map 

U(x) : M 3 -> SU(N) (3c) 
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M 3 being the compact three-manifold on which the theory is defined. Thus, although the 
action is not gauge invariant, the generating function would be so if the dimensionless 
Cher n- Simons coupling k is quantised to be an integer. 

Obviously, a mass term A^A^ for the gauge field is gauge non- invariant. Let us, 
however, consider the following term 

C m = -mTr[(A^ + iV~ 1 d lt V)(A' i + iV^d^V)] (4) 

where the SU(N)-valued auxiliary field V is designed to transform in such a way that 
the whole term is gauge invariant, m has dimensions of mass and the auxiliary field 
V is dimensionless. The total action S = J d s x C = J d 3 x (Ccs + Cm) is then gauge 
invariant under the simultaneous gauge transformations 

A M -> A% = U-^AfJJ - iU-^U (5a) 

and 

V — > V u = VU (56) 

In S only terms linear and quadratic in Aq appear. We may therefore readily do the 
A Q integration in the path integral after doing a Wick rotation t — > it. After the A 
integration, the lagrangian density takes the form 



C = -- [{A* + (iV-'diVrXAa + (iV-'d^a) + (V-'doVnV-'doVU 
k2 f ™ , 45rm /v-i fl m«i[ P , 47rm /T^-ifl m i k 



m + (V-'doVrmja + (V-'doVU - —e*>A?6 ab doA 



b 



327r 2 m L 13 ' k v u ; JL lja ' fe v u ;ai 8tt 
The in the above equation is the usual field strength tensor defined by 

F% = - d s A$ - flA\A) 

In order to be able to pick the relevant degrees of freedom of each of the fields appearing 
in the above terms in the infra-red limit, it is useful to put the system in a box. The 



dimensions of the box L will subsequently be taken to infinity. Introduce, therefore, 

... 1 ^-^ 27rix ■ n 
Ai (x, t) = j exp q in (t) (6a) 

and 



L 

n 



V(x,t) = ^>xp — v n (t) (66) 

n 

In the long wavelength, or the infra-red, limit all the massive modes will be suppressed 
and we can approximate the sums on the right hand sides of the above equations by 
just the zero momentum components. In what follows, the spatially constant modes 
thus obtained for A^ and V are simply denoted by and v respectively. If we now take 
the limit L — > oo in the action and drop all the irrelevant terms we get 

S = Jdt \-^q*V ah q) - ^q?q ia - m^ 1 d v) a (v- 1 d v) a ] (7a) 

where V is the covariant derivative in the adjoint representation 

V ah = 8abd Q -fZ b {v- 1 d Q v) c (76) 

Introducing local coordinates £ a on the group manifold we now have 

Tr(v- 1 d v) 2 = 9ab (0i a i b 

where g a b(Q is the metric induced on the manifold by the map in Eq.(3c), and is given 
by 

9^)=Tr(v-'^-^) (8) 
The action now takes the form 

S = Jdt \-^q-V ah q) - ^q?q ia - 2m^(0C^] (9a) 

The covariant derivative in the above equation is written in terms of £ a as 

V ab = 5 ab d - f c ab g cd (m d (96) 



Let us first consider the ^-independent part of the Lagrangian which we denote by L . 

L o = ~^ ij< ltQja ~ ^Qfqia (10) 

The momentum conjugate to q is then given by 

dL k 

p > = m- = s £y,, ° <lla) 

The above equation yields the following second class constraints 



It is straightforward to evaluate the Dirac bracket [9] for such a constrained system and 
it yields 



47r rab 



{q?,q b j }DB = T e ji 5 ab (12a) 



We can now work out the canonical momentum as 



An identical expression is arrived at directly by using the Faddeev-Jackiw formalism 
[9] . The Hamiltonian is given by the usual Legendre prescription as follows 

Htt= M-L=«%)>(m H *.y«£) (13) 

if we eliminate one of the degrees of freedom using the definition of the canonical mo- 
mentum in Eq.(12b) to arrive at the second equality. Thus, the ^-independent part of 
the Hamiltonian just represents a bunch of (N 2 — 1) harmonic oscillators. This is very 
much like in abelian Chern-Simons quantum mechanics of [5,6]. It is curious, though, 
to note in passing, that the non-abelian feature of the problem is reflected in the fact 
that the frequency of each of these oscillators is quantised in units of 
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We now turn our attention to the more interesting ^-dependent part of the La- 
grangian which we denote by L^. 

Ls = -2mg ab (0i a i b + ^ ^VaWC)^ (14) 
Once again, the momentum canonically conjugate to £ can be written as 

t° = ^ = -4m? + ^Q^fbc (15) 
The corresponding Hamiltonian works out to be 

H * = ~ : L 9ab{0va ' pb (i6a) 

where 

V a = \{^-^ q t q d 3 f a cd ) (166) 

Physically, this part of the Hamiltonian just represents the dynamics of a particle moving 
on the SU(N) group manifold. The extra term in V a is, however, reminiscent of the 
gauge field in the magnetic translation operator of the Landau problem [8]. Indeed, 
in this case it leads to a gauge invariant coupling with the non-abelian Chern-Simons 
particle whose Hamiltonian is given in Eq.(13). It may be appropriate to mention at 
this stage that H% commutes with H . Since the latter is trivially solved, we shall 
henceforth concentrate our efforts exclusively on the ^-dependent part of the theory. 

To get a better insight into the problem, it is prudent, at this juncture, for us to 
specialise to the case of SU(2). An element v G SU(2) can always be represented in 
terms of Euler angles a, 0, 7 and the Pauli matrices as follows: 

CK f3 7 

v = exp(«x 3 -) exp(«r 2 -) exp(«r 3 -) (17a) 

where 

< a < 2tt, < /3 < 7T, -2tt < 7 < 2tt (176) 
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9ab = - J I 1 I (19) 



Now let 

e = a , e = p, e = 7 (is) 

The metric on the SU(2) group manifold can then be calculated using the formula in 
Eq.(8) and works out to be 

1 I 1 cos/3 
' ! 1 

cos/3 1 

The determinant of the metric is easily shown to be 

det | g ab |= -^sin 2 /3 (20) 

The metric is therefore singular at (3 = and (3 = it. These two points, for which the 
analysis needs to be done separately, will be ignored in the following. Substituting the 
metric in the expression for we get 

Hz = ^-{VabV a V b + cos p^V 3 +V 3 V 1 )} (21) 

where rj ab is the flat metric. Using the fact that for SU{2), f abc = e abc we may recast 
the expression for V as 

V a = -(iv a -2L a ) (22) 

where L a are the components of the angular momentum operator. It is useful to record 
the commutation relations satisfied by the operators in 

[e, e 6 ] = k a , A = o, [e, A = in g a \o (23a) 

and 

jf, 

[V a ,V b ] = ihe abc L c = —e abc (n c - 2P C ) (236) 

It is interesting to note that the system has three degrees of freedom viz., a, 0, and 7. 
However, a and 7 are cyclic coordinates. Hence, it also has three constants of motion, 



namely, 7r a ,7r 7 , and the Hamiltonian, H^. It is therefore integrable. In fact, since V a 
depends on the angular momentum operator, L a , we have one integrable system for 
each representation of the Casimir L 2 . This is the infinite class of quantum integrable 
models advertised in the title of the paper and alluded to earlier. 

Since 7r a and 7r 7 are constants of motion, let us, without further ado, set them 
equal to zero. The Hamiltonian, henceforth to be simply denoted by H, reduces to 

H = 4^ - h L > + iL^ + 4^< L ^ ( 24 » 

where we have identified L 1 , L 2 , and L 3 with the components of the angular momentum 
L y ,L 2 , and L x respectively and dropped the suffix on the second component of the 
momentum for elegance. With L 2 and L z being diagonal, as usual, we get for the 
Hamiltonian in the singlet representation, H s , 

H - = ^ (25 » 

i.e., it is just the Hamiltonian for a free particle. The triplet representation presents the 
simplest non-trivial model in the infinite class of integrable models presented above and 
will now be discussed. Introduce then the angular momentum matrices for the triplet 
representation of SU(2) group. 



Plugging in these matrices in the equation for H we get for the Hamiltonian in the 
triplet representation, H t , 

(mR-*R -I^os/^ 

it | (20) 

\l^os/3 ^ + ^ 
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The Schrodinger equation can then be written as 







The equation for ip 2 



is easily solved to get 





16m M /3 





-ik cos P \ 



h 2 ^2 _ ih 2 



16m /3 4m 




16m 



4^2 (P) = eMP) 




4(3 4/3 

ip2 = Ci exp(i— vme) + C2 exp(— i— yme) 

I L lb 



(27) 



(28a) 



(286) 



C\ and C 2 being arbitrary constants. This is very much like in the Landau problem [8] 
where the component of the wavefunction along the direction of the magnetic field is 
just given by a plane wave solution. The non-trivial physics resides in the plane perpen- 
dicular to the magnetic field. It is also well-known that, in this plane, the problem can 
be mapped into the problem of a particle in a harmonic oscillator potential. As already 
mentioned, the present problem is a non-linear generalisation of the above situation in 
the sense that the harmonic oscillator gives way to the mathematical pendulum. Fur- 
ther, because of the non-abelian nature of the problem, we have a matrix mathematical 
pendulum to contend with. With this understanding we may now return to the case at 
hand i.e., obtaining the other two components of ip((3) by solving the two simultaneous 
differential equations 



where 



(4 -4za»Vi(/3) +4icos/?V 3 (/?) = E^{(3) 



E 



16m 



(29a) 
(296) 

(30) 
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This is a set of second order linear differential equations with periodic coefficients. Hence 
the solution we are looking for can be expanded in a Fourier series as 

oo 

ipi (0) = a n exp(m/3) (31a) 

n= — oo 
oo 

MP)= KeMinp) (316) 

n= — oo 

where the coefficients a n and b n are, in general, complex. It is worth mentioning here 
that the solutions of a system of differential equations with periodic coefficients need 
not, in general, be periodic [10]. Indeed, for a generic period, the wavefunction is only 
required to return to itself modulo a phase which is determined by an appropriate 
unitary representation of the translation group. By setting this phase equal to the 
identity we are effectively restricting ourselves to the centre of each energy band in the 
Bloch picture. It is only for this subset of eigenstates that we will be able to obtain 
some exact results. Substituting the above expansions in Eq.(29) we get, after collecting 
terms, and setting the coefficient of each Fourier mode to zero, 

-f(n)a n + 2i(b n _ 1 +b n+1 ) =0 (32a) 

and 

-f(-n)b n - 2z(a n _i + a n+1 ) = (326) 

where 

f(n) = n 2 - An + E (32c) 

The above two equations can be recast as a single matrix equation in terms of a two 
component column vector A n . 

*• = (£) (33) 

After some algebra, this equation can be written in the following compact form 

A n _i + A n+l + M n A n = (34) 
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the matrix M n being given by the equation 

M f(n) + f(-n) f(n)-f(-n) 

M n = cr 2 H a 2 a 3 (35) 

Before we embark on solving this equation, it is worth recalling what we are interested 
in. In terms of A n , the two components of the wavefunction we are trying to solve for 
can be assembled into a column vector, which, for convenience, we shall denote by ip. 
We therefore have 

oo 

ij)(P)= A n ex-p(in(3) (36a) 

n= — oo 

and 

oo 

^((3)= 4>exp(-in/3) (366) 

n= — oo 

Notice now that the system of Eq.(29) is invariant under the discrete transformations 

/?--/?, Vl-V>3, ^s^-^i (37) 

Hence, if ipi(/3) and ^((3) is a set of solutions of the above equations, so is the set 
ip3(—/3) and —ipi(—(3). Once again, as in the discussion following Eq. (31), the latter set 
could differ from the former by a constant multiplicative matrix. If we restrict ourselves 
to solutions which are strictly symmetric under the above discrete transformations, 
however, we have, 

This can be further translated into a relation between A n and A— n as follows. 

A_ n = -ia 2 A n (386) 
Substituting the above results in Eqs.(36) we get 

oo 

i>(0) = - i&2)A n cosn/? + i(l + ioi)A n smn/3] 

n=0 
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and 

oo 

It easily follows from the above two equations that 

„ n oo oo 

/ V t W(/5)^ = 47rJ]^A n = 47rJ](|a n | 2 + |6 n | 2 ) (39) 

^° n=0 n=0 

We shall argue that the above sum is, in general, convergent. Hence the wavefunc- 
tion is square integrable. We shall also derive closed expressions for the above set of 
eigenf unctions. 

In Eq. (32a, b) let us eliminate one variable, say b n . The resulting equation 
for a n couples a n with a n -2 and a n+ 2- This is significant because the even and odd 
coefficients completely decouple and may be treated separately. Since we are trying to 
solve a three-term recursion relation we need to fix two coefficients. We may set all the 
odd coefficients to zero by choosing the first two of them to vanish. Let us then deal 
only with the even coefficients. It is straightforward to show that these coefficients obey 
the recursion relation 

+ /"pap+i + Ppa p = (40a) 



where 



and 



^ = /(-(2 P+ D) {m) 



J(-(2p + 1)) + /(-(2p - 1)) - f(2p)f(-(2p + l))/(-(2p - 1)) 
PP 4/(-(2p-l))/(-(2p+l)) (4UCJ 

and 

a p = a 2p (40d) 
If, on the contrary, we choose to set all the even coefficients to zero, the odd coefficients 
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satisfy a similar equation with slightly modified parameters [i and p. Identical consid- 
erations apply for the coefficients b n . 

We now introduce a one (complex) parameter family of the above three-term recur- 
sion relation as shown below and appeal to two theorems * in the theory of orthogonal 
polynomials. The proofs of these theorems, well-known in mathematics literature, can 
be found in [11]. 

Theorem 1: For any non-real A, there exists at least one solution {a n (A)}g° of the 
equation 

(A) + ^ n a n+ i(A) + p n a n (X) = Xa n (X) 

for which 

oo 

^ I «n(A) | 2 < OO 
n=0 

The above equation actually has two linearly independent solutions by virtue of 
the fact that it is a second order linear difference equation. These solutions P n (A) and 
Q n (X) which satisfy the conditions 

P (A) = 1, P 1 (A) = ^^ 

and 

Q o (A) = 0, Qi(A) = — 

are called orthogonal polynomials of the first and second kind respectively. Hence, let 
a n (X) be given by P n (X). Then, 

Theorem 2: If the series J2^Lo I Pn(X) | 2 converges at any no n- real point A, then it 
converges uniformly in every finite part of the complex A-plane. An identical result also 
holds for Qn(X). 



* These theorems, due to E. Hellinger and R. Nevanlinna, are discrete analogues of 
Weyl's theorems concerning Sturm-Liouville differential equations on the semi-infinite 
real line. 
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It then follows immediately, by choosing the finite part of the complex plane to be 
the origin, and a n to be one or the other of the two types of orthogonal polynomials 
mentioned above, that Y^=o I a n(0) | 2 < oo. The convergence of the series Yl I Ai(0) | 2 
can be argued analogously. Note further that for A = 0, the one parameter family of 
recursion relations we have introduced collapses to the recursion relation we are trying to 
solve. Since identical considerations also apply for 6 n s, we may conclude that the sums 
on the right hand side of Eq. (39) are in general convergent and hence, the wavefunction 
is square integrable. 

The a n or (3 n , hence a n , and similarly b n , and consequently A n , can be obtained 
in a closed form. In order to do this let us consider a related problem. This problem is 
specified by an infinite sequence of orthonormal orbitals {uo, u\, 112, • • •} and a set of real 
parameters {po? Pi-, Pi, • • • /-*o> A*i, 1^2, • • •} which describe the action of the Hamiltonian 
Tt on the orbitals by a symmetric three-term recursion relation [12, 13] 

Tiu n = p n U n + HnU n+ l + Hn-xUn-X (41) 

The a n will be related to the eigenstates of the above Hamiltonian in a very specific 
way. The chain model defined above has a nice physical interpretation. The orbital uq 
represents the initial state of the system and it could be, for example, an electron on 
a particular atom in a solid. The chain model Hamiltonian allows not only a certain 
on-site probability for each of the orbitals that are available, but also an amplitude for 
hopping from one site on the one dimensional lattice to the nearest and the next nearest 
neighbour sites. The example we are trying to solve is a particularly complicated chain 
model in the sense that both the hopping amplitude and the on-site amplitude are 
inhomogeneous. Let us, however, press ahead with our programme, undeterred, for a 
while. We begin by introducing a basis for the orbitals u n . This basis is specified by 
representing u n as an infinite component vector all of whose elements except the nth 
are zero. The non-zero element is, further, chosen to be one. Remembering that n runs 
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from to oo, then 



u = 







W 



Ui 



1 



W 



(42) 



It follows rather easily that the Hamiltonian H is expressed in this basis as a Jacobi (or 
a tridiagonal) matrix 



n 



/ po Po 
Po Pi Pi 

Pi P2 Pi 



\ 



(43) 



V • •/ 

In this representation it is obvious why the chain Hamiltonian is called a symmetric 
chain Hamiltonian. An eigenstate of the Hamiltonian is some linear combination of the 
states {uq, u\, • • •} denoted below by a deliberate misuse of notation as 



n=0 



such that 



H(f> = \<p 



(44) 



(45) 



The above time-independent Schrodinger equation can be written in a matrix form as 

/(A-po) -Po \ 
-Po (A -pi) -pi 

-^i (A-p 2 ) -P2 
-P2 (A - p 3 ) 



a 2 



= 



(46) 



V • / 



V •/ 

The eigenvalues of the Hamiltonian H are given by the zeroes of the determinant A (A) 
of the above matrix. The Jacobi form of the above matrix gives rise to a very simple 
recursion relation for its determinant. It is easy to show that if we define A_i(A) = 
and A (A) = 1 then, 



/^nA n -|-i(A) + p n _iA n _i(A) + p n A n (A) = AA n (A) 
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(47) 



The determinants appearing in the above equation are polynomials in A with the same 
order as their subscripts. If we now substitute for in Eq.(45) and collect coefficients 
of each orthogonal orbital u n , we get 

At n _iex£i + Mn«i+i + Pna^ = X p a^ (48) 

where X p is a zero of Ajv(A) and hence an eigenvalue of the iV-chain, N being a generic 
integer which will be finally taken to oo. affi defines an eigenstate corresponding to an 
eigenvalue X p through the relation 

By comparing Eq.(47) and (48) we get, apart from a common normalization, 

a? 1 = A n (X p ) (49) 

If we now assume that there always exists an eigenvalue X p = 0, 

a n = A n (0) (50) 

yield solutions of Eq.(40a) that we had initially embarked upon solving. A few words 
regarding the properties of the above determinants are now in order. First, it is obvious 
that the secular determinants of the type discussed above would be finite for any model 
with iV points on the chain where N is a generic, but finite integer. Since we are finally 
interested in solving an infinite chain model, we need to ensure that the determinant of 
the infinite dimensional matrix under consideration converges. This is, in general, not 
possible. However, if one plots the zeroes of successively larger determinants, it is known 
that in the infinite limit, they converge to the eigenvalues of the infinite chain. That 
is to say that the determinant converges on the spectrum of the theory but diverges 
for A lying between two energy levels. For getting the /3 n s one has to repeat the entire 
exercise above with the parameters p and p appropriately modified. 
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We shall now, as promised, sketch a method for determining the eigenvalues of the 
above model. Towards such an end, let us recast Eq. (34) as 

*r)= *■*(£.) (6ia) 

where R n is a block upper triangular matrix defined as follows: 

^(-y E) - 1 ) («») 

It follows by iterating the map in the above equation that 

A n+1 \ _ ^ (A x 



where 

T n = R n R n _i ■ • • Ri (526) 

This is as it should be. The above equation for A n merely states that we can work 
out all the coefficients if two of them are chosen and if the 4x4 transfer matrix T n 
is known. The matrix T n , like the matrices R n , R n _i, • • Ri, has elements which are 
functions of the variable E, which could, at random, take any value in the spectrum of 
the theory. It is well-known that this random matrix satisfies the Oseledec condition 
[14] given by 

lim (Tj/Tj^ = S _1 DS (53a) 

where 

D = diag (exp — Ai, exp — A 2 , exp — A 3 , exp — A4) (536) 

The Aj are called the Lyapunov exponents of the infinite chain model under consid- 
eration and are functions of the random variable E. If we choose periodic boundary 
conditions, the matrix on the left hand side of Eq.(53a) is one of the many 2nth roots 
of the identity matrix i.e., it is described by a set of unitary matrices. As n — > 00, this 
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set of matrices densely fills up a limiting unit circle. It then easily follows that S is also 
a unitary matrix and 

[lim(T^T n )*]t = S- 1 D*S (53c) 

n^oo 

Taking the product of the matrices on the left hand sides of Eq. (53a) and (53c) and 
equating it to the product of the matrices on the right hand sides we get, 

DD* = 1 (54) 

It follows, therefore, that the elements of the diagonal matrix D are pure phases. This 
information can now be fed back into Eqn. (53a). The determinant of D is purely 
imaginary and is given by the determinant of a unitary matrix on the left hand side of 
Eqn. (53a) . The latter, however, is not unique. If we pick the simplest root of identity 
i.e., the identity matrix itself, on the densely filled limiting circle introduced above, the 
matching of the determinants yields 

3f(Ai + A 2 + A 3 + A 4 ) = 2nr (55) 

where r is an arbitrary integer. The above equation restricts the randomness in E to 
some specific functional form in terms of r which may be identified with the spectrum 
of the theory. For other choices of the root of identity, the right hand side of the above 
equation gets shifted by a real constant. It is difficult to make progress beyond this 
analytically without knowing the Lyapunov exponents. 

We will conclude this paper by summarising the main results and outlining some 
prospective directions for further research. Amongst the principal results of this paper is 
the fact that the infra-red, or the zero-momentum, limit of massive non-abelian Chern- 
Simons theory is described by an infinite class of models whose elements are labelled 
by the angular momentum eigenvalue. The first non-trivial model in this class was 
examined in considerable detail and closed expressions for some of the eigenstates were 
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obtained in terms of some generalised orthogonal polynomials. A method which, in 
principle, allows us to find the spectrum of the theory was also sketched. Physically 
this model corresponds to the coupling of a non-abelian Chern-Simons particle with a 
quantum mechanical non-linear sigma model or, equivalently, to a particle moving on 
the three-sphere - the SU(2) group manifold. It would be interesting to see what the 
models in the infinite class presented by us correspond to physically for higher angular 
momentum eigenvalues. The two singular points on the three-sphere corresponding to 
(3 = and = ir, that we have ignored in our analysis, merit discussion in their own 
right. One has to probably use a more refined mathematical method, where one chooses 
more than one coordinate chart to cover the manifold, in order to study this problem. 
It is also probable that following the lines of [6] one may be able to derive a phase- 
space path integral localisation formula for the above model. What is perhaps most 
fascinating is the connection that our research opens up between Chern-Simons gauge 
theories and Anderson's chain models. These models, which are of perennial interest 
in the study of alloys, quasi-crystals, and other disordered systems, are a subject of 
intense investigations in condensed matter physics [12, 13]. The assumption we made in 
Eq. (50) regarding the existence of a zero eigenvalue for the chain model is important 
enough to deserve further scrutiny. There exist methods for testing whether the density 
of states in such models is peaked around zero, but the evidence in this regard comes 
from not-so-illuminating numerical work. Finally, one can't help wondering if there 
is a relationship between our model and the so-called quasi-exactly soluble models of 
quantum mechanics [15] that have generated some excitement lately. It may be recalled 
that a quasi-exactly soluble quantum mechanical model is one in which the model is 
exactly soluble for a certain range of a continuous parameter on which the Hamiltonian 
depends. In our case, instead of a continuous parameter, the Hamiltonian depends on 
an infinite, but discrete, set of integers which are the angular momentum eigenvalues. 
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Although we have studied a non-trivial model associated with only one value in this 
set, for higher angular momentum eigenvalues, as well as for groups other than 577(2), 
we still expect a discrete spectrum and square integrable wavefunctions because in each 
case the Hamiltonian is described by a self-adjoint operator on a compact manifold. It is 
however not clear how one obtains the eigenstates in these more complicated cases in a 
closed form. Indeed, this might well be impossible. We may therefore conclude that our 
model is yet another example - clearly more complicated than others hitherto known - of 
a quasi-exactly soluble model. Indeed, at hindsight, the commonality of the important 
role played by the mathematics of orthogonal polynomials in our research and that of 
the subject of quasi-exactly soluble quantum mechanical models seems more than a 
mere accident. 
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